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Introduction 


In a previous study, severe acute respiratory syndrome coronavirus (SARS- 
CoV) was cultured in the presence of bananin, an effective adamantane- 
related molecule with antiviral activity. In the present study, we show that 
all bananin-resistant variants exhibit mutations in helicase and membrane 
protein, although no evidence of bananin interference on their mutual 
interaction has been found. A structural analysis on protein sequence 
mutations found in SARS-CoV bananin-resistant variants was performed. 
The S259/L mutation of SARS-CoV helicase is always found in all the 
identified bananin-resistant variants, suggesting a primary role of this 
mutation site for bananin activity. From a structural analysis of SARS- 
CoV predicted helicase structure, $259 is found in a hydrophilic surface 
pocket, far from the enzyme active sites and outside the helicase dimer 
interface. The S/L substitution causes a pocket volume reduction that 
weakens the interaction between bananin and SARS-CoV mutated helicase, 
suggesting a possible mechanism for bananin antiviral activity. 


At the beginning of the 21st Century, a novel virus, 
the severe acute respiratory syndrome coronavirus 
(SARS-CoV), moved into the human population caus- 
ing SARS with a high rate of mortality. Although the 
last reported epidemic of SARS dates back to April 
2004, the fact that this virus can replicate in a large 
number of animals, including dogs, cats, pigs, mice, 
ferrets, foxes, monkeys and rats [1-3], in addition to 
the natural hosts comprising Chinese palm civets, 
raccoon-dogs and bats [4-6], is of particular concern, 
suggesting that preparedness with vaccines and anti- 


Abbreviations 


viral drugs against this potentially re-emerging agent is 
necessary. 

It has been shown that treatment with ribavirin and 
corticosteroids as possible drugs against SARS-CoV 
only had slight beneficial effects or even enhanced viral 
replication in mice [7,8]. Thus, development of new 
anti-SARS drug is urgently needed for the potential 
SARS re-emergence. The relative conservation and 
essentialness in functionality of a particular gene are 
used as indicators to evaluate a drug target. On the 
basis of these criteria, helicase is a good target, being a 
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relatively conserved protein in SARS-CoV (e.g. a less 
variable protein compared to spike protein) and criti- 
cal for viral replication [9]. Accordingly, the latter pro- 
tein has been proposed as an attractive target for anti- 
SARS research [10], in analogy with the promising 
results obtained for herpes simplex virus-! [11,12]. 

In the present study, a structurally driven investiga- 
tion for the design of new SARS-CoV helicase inhibi- 
tors is performed by correlating the predicted enzyme 
structure [13] with the observed bananin activity 
against SARS-CoV [14]. Bananins are a class of com- 
pounds with a unique structural signature incorporat- 
ing a trioxa-adamantane moiety covalently bound to a 
pyridoxal derivative to add potential cytoprotective 
functionality [15]. Several parent adamantane deriva- 
tives are already used clinically [16,17], although the 
antiviral activity of the newly developed bananin has 
not yet been investigated extensively. Jn vitro assays 
demonstrated that bananin effectively interferes with 
SARS-CoV ATPase activity by inhibiting helicase 
activity. Furthermore, in a cell culture system of 
SARS-CoV, bananin inhibited viral replication at a 
half maximal effective concentration of less than 
10 4M and a concentration causing 50% of cell death 
of over 300 uM, suggesting that it represents a promis- 
ing anti-viral drug candidate [10]. To investigate bana- 
nin primary targets, bananin-resistant viruses were 
selected by culturing SARS-CoV (GZ50 strain; Gen- 
Bank accession number AY304495) on fetal rhesus 
monkey kidney cell line FRhK-4 in the presence of 
high concentrations of this adamantine derivative. The 
half maximal effective concentration of bananin on 
these mutant viruses was demonstrated to be more 
than 50 uM. Mutations were found in helicase (S259L), 
membrane protein (A68V and R124W) and spike pro- 
teins (N479]). The trans-expression of mutant helicase 
or membrane protein during wild-type virus infection 
can rescue viral replication in the presence of bananin, 
demonstrating that the SARS-CoV helicase and 
M proteins were effective drug targets [14]. 

The present study describes the systematic search for 
those mutations found in bananin-resistant SARS-CoV 
variants. Subsequently, structural and functional 
results are compared to define the possible mechanisms 
of bananin activity, and are also used to drive 
restrained docking simulations of the bananin-SARS- 
CoV helicase interaction, aiming to define the sterical 
requirements of new antiviral drugs. 


Results and Discussion 


All the mutations found in the isolated bananin-resis- 
tant SARS-CoV variants are summarized in Table 1. 
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Table 1. Mutations in bananin-resistant virus variants. S, spike 
protein; M, membrane protein; NT, not tested. 


Variants Helicase S M 
B15 S259L N4791 A68V, R124W 
B18 S259L, L297L/F NT A68V, R124W. 
B6, B14 S259L NT A68V, R124W 
Other nine variants S259L NT NT 


Of primary relevance is the fact that the $259L muta- 
tion in helicase is always present, as well as the A68V 
and R124W mutations in membrane protein. This find- 
ing initially suggested that the observed antiviral activ- 
ity could arise from bananin interference on a 
hypothetical helicase-membrane protein interaction, as 
in the case of the closely-related coronavirus mouse 
hepatitis virus, where the replicase protein complex 
including helicase colocalizes with M in the endoplas- 
mic reticulum-Golgi intermediate compartment for vir- 
ion packaging [18]. Thus, co-immunofluorescense and 
co-immunoprecipitation experiments were performed, 
although no evidence of helicase-membrane protein 
interaction could be obtained (Doc. S1 and Fig. S1). 
Thus, the mechanisms of bananin activity have been 
ascribed to the binding of the small molecule to single 
viral proteins and, in particular, to those exhibiting 
mutations in the bananin resistant SARS-CoV variants. 

As shown in Table 1, mutations in spike protein, heli- 
case and membrane protein have been detected. By 
using SCORECONS software [19], which quantifies residue 
conservation in multiple sequence alignments, Shan- 
non’s entropies have been calculated for each of the 
SARS-CoV protein sequence positions where mutations 
have been observed. In the case of SARS-CoV spike 
protein sequence, Shannon’s entropy (in the range from 
0 for invariant to | for hypervariable protein sequence 
positions) reveals that the 479 position, where the N/I 
mutation is found, corresponds to a highly variable site. 
A Shannon’s entropy value of 0.72 is obtained by 
retrieving all of the 92 complete sequences of SARS- 
CoV spike protein present in the National Center for 
Biotechnology Information (NCBI) databases. The fact 
that bananin does not target on SARS-CoV entry 
(Doc. SI and Fig. $2), suggests that no bananin-related 
activity can be attributed to the N479I mutation. 

For the SARS-CoV membrane protein, only second- 
ary structure predictions can be obtained, limiting a 
detailed structural interpretation of the functional roles 
of A68V and R124W mutations. BLAST analysis on the 
19 complete sequences present in the NCBI databases 
for the SARS-CoV membrane protein suggests that the 
observed conservative substitution A68V is also encoun- 
tered in two native viral clones (i.e. dbj_BAE93405 and 
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gb_AAP33701). From a structural point of view, the 
A68V mutation, occurring in a predicted protein trans- 
membrane region, should exhibit only a minor func- 
tional relevance. The case is different for the RI124W 
mutation, which is outside the trans-membrane seg- 
ments, with the arginyl residue being very conserved in 
all the available related sequences. Thus, the zero Shan- 
non’s entropy value calculated for the totally invariant 
124 position of the membrane protein sequence suggests 
some functional role for such R/W replacement. How- 
ever, the absence of tertiary structure information for 
the SARS-CoV membrane protein prevents further 
functional analysis of the R124W mutation. 

The fact that tertiary and quaternary structures can 
be predicted for SARS-CoV helicase [13] allows a 
deeper insight into the function/structure correlations 
for native and mutated forms of the viral enzyme. 
Preliminary analysis on the frequency of amino acid 
substitutions among the 78 complete sequences of 
SARS-CoV helicase available from the NCBI databases 
indicates that S259 and L297 are totally conserved sites. 
Therefore, it can be suggested that chemical pressure as 
a result of the presence of bananin in cell cultures is the 
only driving force for selecting the observed mutations. 

Functional validation for the helicase tertiary struc- 
ture shown in Fig. | (atomic coordinates are available 
from the Protein Model Databank at http://www.caspur. 
it/PMDB under the accession number PM0076418) is 
provided by the observation that, in studies performed 
in vivo, S259/L and L297/F mutations do not interfere 
with viral metabolism because mutant virus variants 
carrying these two point mutations exhibited normal 
replication, as in the wild-type virus [14]. Of primary rele- 
vance is the fact that both $259 and L297 are predicted 


Fig. 1. Predicted quaternary structure of SARS-CoV helicase dimer. 
On each monomer surface, and colored with different gray scales, 
metal-binding domains, DNA duplex, ATP and bananin have been 
highlighted, respectively, in cyan, green, red and yellow. 
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to be outside the surface regions where DNA binding, 
NTPase activity and dimerization occur. Moreover, it is 
interesting to note that L297, replaced by a phenylalanyl 
residue, is totally buried in the helicase structure, and 
that the conservative L/F substitution [20,21] does not 
cause any major changes in the helicase core structure. 

A. structural comparison of molecular models 
obtained for wild-type SARS-CoV helicase and the 
two bioactive mutants indicates how the replacement 
of S259 with the leucyl bulky side chain determines a 
volume decrease of a hydrophilic pocket present on 
the helicase surface. This surface pocket, formed by 
N257 and 1258 backbone atoms together with S259, 
D260 and E261 side chains in the wild-type helicase 
(Fig. 2), reduces its volume from 253.45 to 211.40 A* 


Fig. 2. Predicted bananin binding pocket of SARS-CoV helicase (A) 
in wild-type and (B) bananin-resistant variants. Seryl and leucyl 
side chains in position 259 are shown in a green ball and stick 
representation. 
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upon S/L replacement. In the case that the latter 
hydrophilic pocket of SARS-CoV helicase is the bana- 
nin binding site, the S/L mutation weakens the inter- 
molecular interaction by reducing the number of 
possible hydrogen bond formations. 

Electrostatic potential analysis for wild-type and 
S259/L mutant forms of SARS-CoV helicase was also 
carried out. As shown in Fig. 3, identical surfaces 
charge distributions are obtained, and therefore no 
electrostatic effects on helicase dimerization or heli- 
case-DNA interaction can be attributed to the latter 
mutation. Furthermore, complementarity of positive 
and negative charges at the dimer interface region is 
apparent, supporting the reliability of the SARS-CoV 
helicase predicted structure. 

Quantitative evaluation of S259/L mutation effects 
on bananin-SARS-CoV helicase binding has been per- 
formed with docking simulations on wild-type and 
mutated forms of the viral enzyme. In Fig. 4, the 
modes of bananin interaction with wild-type helicase 
are shown according to the results obtained from the 
docking simulation procedure. Thus, it is apparent 
how surface pockets of native and S/L mutated heli- 
cases are differently filled by bananin because binding 
with the small molecule involves a larger molecular 
interface in the case of the former helicase. Further- 
more, the absence of the S259 OH group in the 
mutated viral enzyme prevents the formation of one 
hydrogen bond with bananin, accounting for the 
reduced strength of the bananin-helicase interaction. 

Explanations of drug activity are usually provided 
by conformational changes in the targeted protein or 
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by competitive binding at the protein active site. Alter- 
native mechanisms of bananin antiviral activity have 
to be found because its protein target at the 259 
sequence position presents a fully exposed side chain, 
and hence only limited local conformational rearrange- 
ments should result from the S/L mutation. Moreover, 
the fact that this mutation site is very far from the 
active site suggests that S259 is located either in an 
allosteric site of the enzyme or in a critical position for 
the overall protein flexibility. The fact that the $259/L 
helicase mutant is fully active is consistent with the 
above hypothesis proposing that this amino acid sub- 
stitution, which is critical for bananin binding, does 
not interfere in the interaction of the enzyme with its 
natural substrates. The presence of the leucyl side 
chain in the protein mutant appears to cause steric 
hindrance to bananin binding, leaving the traffic of 
water molecules in and out from the hydrophilic 
pocket almost unaffected. Removal of such water mol- 
ecules upon bananin binding could reduce helicase 
flexibility, which a very critical feature for the activity 
of this class of enzymes [22]. 

Thus, it can be concluded that bananin resistant 
SARS-CoV variants have delineated an overall protein 
mutation pattern indicating the critical role of the heli- 
case $259/L mutation. The possibility that bananin 
binding to S259 reduces the enzyme activity affecting 
helicase dynamics is consistent with the observed bana- 
nin-resistance of SARS-CoV variants containing 
S259/L helicase mutations. These results are useful 
with respect to the rational design of new anti-SARS- 
CoV drugs in the event of a new unexpected pandemic. 


Fig. 3. Electrostatic potential distribution 

of SARS-CoV helicase: the basic (dark blue 
coloured) region involved in DNA binding 
and the dimer interface (circled in green) are 
shown. On the opposite side of the protein, 
S$259/L regions in the wild-type (inset 1) 
and mutated form (inset 2) of the viral 
enzyme are also shown. 


386 FEBS Journal 278 (2011) 383-389 © 2010 The Authors Journal compilation © 2010 FEBS 


Z. Wang et al. 


Fig. 4. Lowest energy structure of the wild-type SARS-CoV 
helicase-bananin complex obtained from the docking simulation (A). 
Bananin hydrogen bond network with helicase donor/acceptor 
moieties are also shown (B). 


Experimental procedures 


Generation of bananin-resistant virus 


SARS-CoV strain GZ50 [23] was cultured on FRhK-4 cells. 
This cell line was used to isolate and culture this virus 
strain from the very beginning, and was considered fully 
permissive for viral replication [10,14]. SARS-CoV was cul- 
tured in the presence of 50 tM bananin for four passages 
and then 100 um bananin for an additional four passages. 
The bananin-resistant virus variants were identified by a 
plaque assay in the presence of 100 UM bananin as 
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described previously [24] and further isolated by isolating 
viral plaques. 


Sequencing of mutant virus genome 


Fourteen pairs of primers were designed according to GZ50 
sequence for PCR amplification of the whole genome of 
mutant SARS-CoV. Each agarose gel purified a fragment 
of approximately 2 kb that was used as a template for the 
sequencing reaction. PCR primers and sequencing primers 
are available upon enquiry. 


Protein sequence analysis 


Protein sequences of SARS-CoV helicase (SP_P59641), spike 
(SP_P59594) and membrane (SP_P59596) proteins were 
retrieved from SwissProt Database. Sequence alignments of 
these three proteins with all SARS coronavirus sequences 
were obtained with CLUSTALW, version 1.8 [25], and analyzed 
in terms of sequence variability by using the SCORECONS ser- 
ver [19]. Shannon’s sequence entropies were considered as a 
quantitative measure of residue conservation. 


Molecular modeling 


The predicted structure of SARS-CoV helicase, taken from 
the Protein Data Bank (PDB) with the PDB ID code 2G1F 
[13], was used as the initial reference structure. By using 
GROMACS software [26], ten cycles of simulated annealing of 
500 ps each were carried out to improve side chain packing 
and to remove most of the stereochemical ambiguities pres- 
ent in the selected PDB file. Similarly, the structures of heli- 
case S/L and L/F mutants were obtained, and the atomic 
coordinates of the lowest energy structures of the wild-type 
form of SARS-CoV helicase are available from the Protein 
Model Databank  (http://www.caspur.it/PMDB) under 
accession number PM0076418. By structural homology with 
other helicase dimers and helicase-DNA adducts (PDB ID 
codes IUAA, 3PJR and 2IS1), the interface between 
SARS-CoV helicase and DNA has also been predicted. The 
structure of bananin (i.e. 1-[3-hydroxy-5-(hydroxymethyl)-2- 
methyl-4-pyridinyl]-2,8,9-trioxaadamantane-3,5,7-triol) was 
parameterized by using MOPAC2007 [27]. Volumes of the 
proposed bananin binding pocket of native and mutated 
SARS-CoV helicases were measured using the online tool 
caAsTp (http://cast.engr.uic.edu) [28]. ADAPTIVE POISSON— 
BOLTZMANN SOLVER software was used for evaluating the 
electrostatic properties of SARS-CoV helicase [29]. Figures 
were prepared with PYMOL using PDB2PQR [30,31]. 


Docking simulations 


AUTODOCK 4.0 was used to simulate a flexible docking pro- 
cess for the interaction of bananin with SARS-CoV helicase 
and to analyze their binding modes [32]. The AUTODOCK 
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tool (ADT) was used to optimize ligand and protein by add- 
ing polar hydrogens and loading Kollman united atoms 
charges, as well as to perform docking calculations. A grid 
box with dimensions 40 x 40 x 58 points was constructed 
around the SARS-CoV helicase S259 residue. All bond 
rotations and torsions for bananin were automatically set 
by the ADT routine. The Lamarckian genetic algorithm pro- 
cedure was employed and the docking runs were set to 250 
and a maximum number of 2 500 000 energy evaluations. 
In the docking simulations, all the other parameters were 
set to defaults. The resulting orientations with rmsd 
< 0.5 A were clustered. 
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Doc. S1. Supplementary material. 

Fig. S1. No interactions between SARS-CoV helicase 
and M proteins. 

Fig.S2. Bananin has no effect on HIV/SCV pseudo- 
typed viral entry. 
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